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. Abstract 

We present a next-to-leading evaluation of the resummed coefficient function for the shape function. The 
' results confirm our previous leading-order analysis, namely that the coefficient function is short-distance- 
C^h! dominated, and allow relating the shape function computed with a non-perturbative technique to the physical 
^ • QCD distributions. 

> 
• 1—1 

X 
J3 



On leave of absence from Dipartimento di Fisica, Universita' di Roma I, Piazzale Aldo Moro 2, 00185 Roma, 
Italy. E-mail address: ugo.aglietti@cern.ch. 



CERN-TH/2001-050 
March 2001 



1 Introduction 



In this note we present a next-to-leading evaluation of the coefficient function of the shape function. The 
leading-order analysis was done in [Q. The coefficient function allows relating the shape function — computed 
with a non-perturbative technique — to physical distributions in semi-inclusive heavy-flavour decays. 

In general, we consider the process 

X + (non - QCD partons) (1) 

in the hard limit 

Q > A (2) 

and in the semi-inclusive region^ 

M < Q, (3) 

where M is the invariant mass of X, the hadronic final state; Hh is a hadron containing a heavy quark h and 
the non-QCD partons can be a lepton pair, a vector boson, a photon, etc. The quantity Q is the hard scale 
of the time-like process: Q = 28, where £ is the final hadronic energy and A is the QCD scale. Well-known 
examples of (0) are: 

B^Xs+j (4) 

for a large photon energy and 

B^Xu + l + i' (5) 
for a small hadronic mass or a large electron energj^. 

For the hadron at rest, i.e. with velocity w = (1; 0, 0, 0) , and the jet X flying along the minus direction (— z 
axis), the shape function is defined as |§|-(5|: 

ip{k+) = {HUv)\hlS{k+-iD+)K\Hh{v)). (6) 

The latter is also called structure function of the heavy flavour and represents the probability that the heavy 
quark in the hadron has momenta 

Ph = ruHV + k (7) 

with any transverse and minus component and with given plus component The static field hy (x) is related 
to the Dirac field of the heavy quark h (x) by: 

h (x) = e-™"''-"/^, {x) + O (—] . (8) 

\mHj 

With the shape function, the decay (|l|) is related to its respective quark-level process 

/i -> X + (non - QCD partons) , (9) 

where the heavy quark has the momentum with the distribution (H). Note that the state X, unlike X, does 
not contain the light valence quark(s) in H. The invariant mass m of X is related to the plus virtuality of the 
heavy quark by the relation: 

= -f . (10) 



'^This region is also called threshold region, large- region, radiation-inhibited region and Sudakov region. 
^For the rare decay (W), since Q = I 1 -|- /"^s ) i one can actually set Q = mg. 
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The shape function is a non-perturbative distribution — analogous to parton distribution functions — and it 
describes the shce of the semi-inclusive region in which 



TO^ - AQ. 

In terms of the variables used in QCD and in the effective theory (ET), this is: 



fc+~A, ^-^1-^' (11) 



where 



III' / ^ - 

^-l-Q,, (12) 

and N is the moment index (see later). 

The shape function describes the mass distribution of a jet coming from the decay of a heavy flavour. It 
therefore generalizes to region ( pi] ) the jet function J {Q^, m^) introduced in ref. representing the probability 
that a light parton produced in a hard collision with scale Q evolves inclusively into a jet with mass m. The 
main difference is that in our case the distribution depends on the process as a whole, and not only on the 
fragmentation of the final quark. The shape function depends on the initial hadron state; for example, it is 
different for a B meson and a Af, hyperon. 

The coefficient function C is defined by the relation^ 

V5 (fc+; Q)= f dk'+ C (fc+ - fcV; Q, ^ {k'+'ii) , (13) 



where ^ is the renormalization point or UV cut-off of the operator entering the definition of the shape function 
(see eq. (||)). The coefficient function is then obtained by evaluating in next-to- leading order (NLO) the QCD 
effective form factor ip{k+;Q) and the shape function (p{k+;fi) and taking their ratio (see later). Since C 
is expected to be a short-distance quantity, we compute the QCD distribution and the shape function in 
perturbation theory (PT) for an on-shell heavy quark (fc = 0). This expectation will be verified a posteriori. 

The coefficient function is different, even at the leading level, for the shape function regulated in dimensional 
regularization (DR) and in lattice regularization; it is short-distancc-dominated in both cases [Q. The main 
difference is that, at one loop, C contains a double logarithm of fc+ in DR, while it contains at most a single 
logarithm of in lattice regularization. In view of the applications, we specialize ourselves in the lattice case. 

The paper is organized as follows. In sec. 2 we present a compact derivation of a resummed QCD form factor 
ip{k+\Q) to NLO accuracy. In sec. 3 we perform a similar computation for the shape function (/?(A;+;/u). In 
sec. 4 we evaluate the coefficient function and we discuss the physical implications. Finally, in sec. 5, we draw 
our conclusions. 



2 The QCD distribution 

We are interested in a general QCD distribution f (z) in the threshold region. For clarity's sake, let us consider 
a specific case: the photon spectrum in the decay (^), 



QCD quantities arc denoted by bold-face symbols, the related quantities in the ET by normal symbols. 
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where 

rub 

and Fb is the Born width: 



(15) 



Tb = r \:T \v..v;f\c, i,b) p. (16) 



The MS mass is related to the pole-mass by: 

rub^rub) = (l - ^) (17) 

The distribution reads, to order as [§: 

i{x) = i{x)+askS{l-x) + asd{x). (18) 
The function f (x) factorizes (perturbative) long-distance effects occurring in (^) and reads:^ 



where 



A, = % = (20) 

TT 4 TT 



and Cf = {N^ ~ l) / (2-?Vc) = 4/3. The spike term involves the constant 



fc = -^flog^ + ? + ^), (21) 



TT \ nib 4 3 

where fib is the renormalization scale of the operator O7; this (unphysical) scale must be taken of O (nib) in 
order to avoid large (ultraviolet) logarithms in the matrix elements. The function d (x) is the "remainder" : 

d{x) = — \7 + x-2x^ -2{l + x) log (1 - x)] . (22) 

Plus-distributions are defined as usual as P{x)^ = P [x) — (5(1 — x) dy P (y) . For the resummation, it is 
natural to write the function f (x) in an "unintegrated" form as 



{{x) = S{l-x)+ de dt 



Aias , Sias , Cias 



et e t 

where we have defined the unitary energy and angular variables 



[(5(l-x-ei)-(5(l-x)], (23) 



E , I-COS6I 

and t^^^. (24) 

The quantity E is two times the energy of the soft gluon, E — 2Eg, and 6 is the gluon emission angle. It proves 
to be convenient to separate the soft from the coUinear term, as we have done. By explicit evaluation, we find 

^1 ^ = (25) 

n An 



Note that Bi = 5*1 + Ci, i.e. it is the sum of the coefficients of the soft and the coUinear terms. Expression (|23| 
is symmetric for e ^ t. 

Next-to-leading corrections are included following these prescriptions: 



^Overlinod quantities denote subtracted quantities, containing only infrared logarithms. 
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1. Replacement of the bare coupling by the two-loop running coupling: 

/ 2^ _ 1 _ Pi logloggVA^ 

"^■^'^^-/3„log,VA2 Pi log^gVA^ ' ^''^ 

evalutated at the gluon transverse momentum squared |^]: 

as ^ as (fci) , (27) 

where 

kl ~ Q^e^t. (28) 

The first two coefficients of the /3-function are: 

nCA-2nF _ 33-2nF 
^" " 12^^ " f2^ ' 
^ _ UCl ~ BCauf -SCpriF _ 153 -19 np 
- 24^^ 247r2 ' 

where Ca — Nc = 3 and n^? = 3 is the number of active quark flavours; 

2. Inclusion of the two-loop correction to the term Aias, so that 

Aias ^ Aias + A2al. (30) 

The explicit computation gives |^ 

A2 = l^K (31) 
where, in the MS scheme for the coupling constant, 

with Tr ==1/2. 
The "effective" one-gluon distribution^ then reads: 

'Aias (Q^H) + A2 al {Q\H) as (QhH) C, as (Q'eH) 



i{x) = 6{l-x)+ [ de [ dt 
Jo Jo 



et € t 

[S{l-x-et)-5{l-x)]. (33) 



This integral is no longer symmetric for e ^ t because of running coupling effects. 

A consistent resummation to NLO accuracy (and beyond) is naturally performed in moment space, so that 
we consider [O: 



IN 



1 _ 

dx x^-^ f [x) 







1 + Afjv. (34) 



^We call it "effective" because the insertion of the running coupling in the time-like region and of the term proportional to A2 
already includes some multiple-gluon-cmission effects. 



With this definition of the Mellin transform, the total rate is given by the first moment. Note that iiy^i — 1. 
The inverse transform to the original distribution in x-space is usually done numerically [l2| . The moments of 
the one-loop distribution f (x) given in eq. ( |l9[ ) read: 



f „ - 1 - 



Aias 



log n — Bias\ogn 
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as + 0{- 



(35) 



We have defined n = N/Nq, where Nq = e'^E = 0.561459 . . . with 7^ = 0.577216 ... the Euler constant g. 
The single logarithm, the logn term, is positive and tends to increase the rate, as is often the case [|l^. Its 
coefficient Bi is rather big; it is over a factor of 2 larger than in DIS, for example, where Ci is the same and 
5*1 = 0. Simple exponentiation of the effective one-gluon distribution takes place in A^-space, so that 



(36) 



The non-logarithmic integrations of the terms proportional to 6*1 and Ci can be explicitly done and one 
obtains: 



log f AT = 



„Af-l 



dx ■ 



- 1 



1-x 
+Si 



1 



dQ r . ( 2^ 

Q2(i-i:)2 q 



Using the large- approximation x 



as (g' (1 - xf) + Cias {Q\1 - x))] + O (ag log""! N) 
(1 — a; — 1/n) |^|^, the moments read: 



N 



(37) 



' s s s s 

logf„ = -TTTT log s log log s - 2 log - log log - -t- log ^ log log ^ 



s s 
log log s - 2 log log - + log log 



log log s - 2 log log - + log log — 



log log s - log log 



log log s — log log — 
n 



(38) 



where s is the square of the hard scale in unit of the QCD scale, 



(39) 



Note that f„ = 1 for n = 1. Substituting logs in eq. (38) by its expression in terms of the two- loop coupling, 

(40) 



log s = -J— + ^ log iPoO^s) , 



one obtains the exponent (|8|) as a series of functions ||i4| : 

f„ exp [L gi iPoasL) + g2 (PoasL) + as gs iPoasL) 
where the leading and ncxt-to- leading functions are respectively [|l5[: 



(41) 



91 (w) 

92 (w) 



A, 1 



2/3o w 
P0A2 - PiA 



[{1 - 2w) log (1 - 2u)) - 2 (1 - w) log (1 - w)] ; 



i [log(l - 2w) - 21og(l - w)] - ^ [log2(l - 2w) - 21og2(l 



2/3g 



+ ;^log(l-2u;)-f ^log(l-H- 



(42) 



®The error in this approximation amounts to next-to-next-to-lcading logarithms and to terms suppressed by powers of l/N. 



5 



We have defined as = as (Q^) ^nd L = logn. The two-loop terms coining from the expansion of the resummed 
distribution read: 



The effective form factor is real — and therefore sensible — only for 

1 



n < Tlcrrt = exp 



(43) 



(44) 



This restriction is induced by the soft terms — i.e. by those proportional to Ai, Si and ^2- The reliability of 
the resummed result is then limited to 



n < 



g 

A 



(45) 



The condition is the same as that in the leading order analysis. On the other hand, the hard coUinear contribution 
— the term proportional to Ci — remains well-defined up to much larger values of n : 



n « "Lt ~ ^ ^ » A2. (46) 

The collinear factor is therefore well-defined in the kinematic region (jl^) described by the shape function. It 
becomes singular only in the resonance or exclusive region rri^ ^ A^, which is well beyond the reach of the shape 
function. 



0.8 - 



0.6- 



0.4 - 




Figure 1: Plot of the QCD form factor f„ as a function of n (first 24 moments) for the values of the parameters 
discussed in the text. Solid line: NLO distribution; dashed line: LO distribution; dotted line: expansion to 
order a| of the exponent in the NLO distribution. 

The LO and NLO form factors are plotted in fig. 1 for Q = ms — 5.2 GeV, for which as ~ 0.21 and Ucrit — 28. 
The NLO curve lies above the LO one and it exceeds one in the small-n region because of the single logarithm, 
as already discussed. The two-loop curve lies above the NLO one because the leading terms in higher orders 
are all negative (cf. eqs. (^) and (ji^)). 

The resummed distribution is usually written as ]l6| : 

iN (as) = K (as) (n (as) + (as) , (47) 
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where K is a, function with a power series expansion in : 

K{as)^l + kas + k'al ~\ (48) 

and D^r is a remainder function: 

DNias)= dNas+d'j^al + --- , (49) 
starting at order as and vanishing in the large- limit, 

Dn- ^ for N ^ oo. (50) 

To NLO accuracy, the contributions of order as to the function K (and Dn) are needed. To the same accuracy, 
in place of (^7|), one can write f_/v — {I + as k + as d n) ^ n ■ 

An expression analogous to ( p7| ) also holds for the distributions in the semileptonic decay (js]) with a fixed 
hadronic energy, i.e. to differential distributions in the hadronic energy. One simply has to make the replacement 

x^z (51) 



in eq. (47) and to account for a dependence of K and Dn on Q/tub- The function f is the same in the 
semileptonic and in the rare decay: 

f(.x)=f(z). (52) 

The general variable for the threshold region of heavy flavour decays is z, even though, for the rare decay (^, 
it is more convenient to use x. In the semi- inclusive region these two variables concide, because 



l-x 

\-z = 



{2-x) 



2 = l-x + o((l-a;)^) . (53) 



3 The shape function 



Let us now consider the quantity in the effective theory related to f (z\ namely the shape function (fc_|_) . The 
latter has the analogous decomposition: 



(fc+) = (fc+) ^- ashb (fc+) + asp (fc+) . 
The function Tpikj^) factorizes long-distance effects and reads: 

-n ^ w, ^ . 0(O;fc+;-/i), -fc+ ^ ^ f" , -Z+ 

^ fc+ b{k+)-Aias ' 7 log i+^ias(5 /c+ / — ^ log ^ 

g(0;fc+;-M) „ , 
+bxas ; Siaso{k+) — — , 

-k+ -1+ 



(54) 



(55) 



where ^ = 2As and we have defined 6 (ai; a2; • ■ • a„) = 9 (ai — 02) (02 — 03) ■ ■ - 9 (a„_i — a„) . The regulariza- 
tion used Q imposes a cut-off on the spatial loop momenta and not on the energies, 

|T| < As, -00 <lo < +00, (56) 
and it is qualitatively similar to lattice regularizatioii[|]. Proceeding in a similar way, one can write: 



dE 

ipik+) = Sik+)+ — dt 



Mas {EH) ^ Mai {EH) ^^^^ 



[5 {k+ + Et) - 5 {k+)] . (57) 



^By this we mean that the double logarithm at one loop is the same in the two regularizations. 
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A few comments are in order. The hard scale Q does not appear in Tp (fc_|_), as it should. The term proportional 
to Ci is absent because the shape function retains only the leading terms in the soft limit i? — > 0, which are 
proportional to \/E. With the same definition of the coupling constant in QCD and in the ET, the constant 
A2 is the same in the two theories jlj, Since A;+ = — Q (1 — z) , the shape function in the "QCD variable" 

z reads: 



f{z) = 5{l-z)+ - dt 



l/r 



de 



t 



Sias {Q^e^t 



[S{l-z-et)-5{l-z)], 
(58) 



where we have defined an adimensional function as 

7{z) = QTp{k+). 

The quantity r is the ratio of the hard scale to the UV cut-off of the effective theory, 

r — — > 1, 



(59) 



(60) 



because the shape function is defined in a low-energy effective field theory. To avoid substantial finite cut-off 
effects, one has also to assume 

/x>A. (61) 
Taking moments and exponentiating the one-loop distribution as in the QCD case, one obtains: 



/ 



N 



(62) 




Figure 2: Plot of the shape function /„ (first 25 moments). Solid line: NLO distribution; dashed line: LO 
distribution; dotted line: expansion to order Og of the exponent in the NLO distribution. 

Making the same approximations as in the QCD case, the shape function is written: 

log 7„ = -9{n- r) f'^ ^1 f'^^^ [A,as {Q'ey) + A^ a| (Qhy)] + S,as (QV) ] • (63) 

Jl/n y [Jy £ J 
®A difference would imply a breakdown of factorization. 
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A straightforward integration then gives: 



log fr, 



{n — r) 



+ 



2/3o 



s s s s s s 

log — log log — - 2 log — log log — + log ^ log log 



+ 



S S S ' 

log log -2 - 2 log log h log log ^ 



+ 



log log ^ - 2 log log h log log 



,s s 
log log ^ - log log ^ 



(64) 




Figure 3: Comparison of QCD form factor with the shape function. Solid line: NLO QCD distribution; big 
dashed line: NLO shape function; small dashed line: LO QCD distribution; dotted line: LO shape function. 

The main point is that the terms containing arc independent of r : the infrared behaviour is the same as in 
QCD, as expected on physical grounds. Substituting log s by its expression in terms of the two-loop coupling, 
the result reads: 



log /„ = L 5f ^ {PqCxsL- PoasR) + {PoasL; p^asR) + 



where: 



gf^ {w; t) = --4^ [(1 - 2w) log (1 - 2w) + (1 - 2t) log (1 - 2t) - 2 (1 - w; - r) log {1-w-t)\, 



9^ {y^; r) = ^3^'^^ [log(l - 2w) + log (1 - 2r) - 2 log {I - w - r))] 



and 



[log2(l - 2w) + log2(l - 2r) - 2 log^ {1 - w - t)] + 
+ ^ [log(l - 2w) - log(l - 2t)] , (65) 



R = logr>0, (66) 
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and the over-all factor 9 {L — R) has been omitted. Note that /„ is continuous at n = r, but it has a cusp 
singularity in this point. The range of n is restricted by eq. (^5|): the same limitation of the QCD distribution 
occurs. The QCD form factor f„, except for the Ci term, is obtained by taking t = in the above expression, 
i.e. letting the gluon energy to reach the hard scale. The shape function is plotted in fig. 2 for the same values of 
the parameters of the QCD distribution and for fi = 2 GeV. In fig. 3 we compare the QCD form factor with the 
shape function. In NLO, the shape-function curve is below the QCD one for small n because it has no Ci term. 
NLO distributions are very close to each other for large n, but this seems accidental; indeed, this phenomenon 
does not occur for the LO distributions. In general, the shape-function curves look somehow shifted with respect 
to the corresponding QCD ones. 

The moments of the shape function can be decomposed in a form similar to the QCD distribution: 

/at (as) = H {as) 7n («s) + Pn (as) , (67) 

where the functions H (as) and Pn (05) have a power series expansion in as : 

H{as) = l + has + h'al + --- 

Pn (as) = PNas+PN^s^ • (68) 

To NLO accuracy, an expression alternative to ( |6^ ) including higher- twist effects is /jv = {I + as h + as pn) In- 



4 The coefficient function 



We introduce the coefficient function Cn relating the shape function to the QCD form factor, as 

fjv ^ Cn In- (69) 

Inserting the expressions previously obtained and neglecting the remainder functions Dn and Pn, the coefficient 
function reads: 

Cn ias)= G{as)^Nias), (70) 

where 

G{as)^^^^^l+ gas + --- (71) 
H (as) 

and 

fjv {as) I 
Sjv {as) = =—. — - (72) 
Jn ("s) 

with g — k ~ h. To NLO accuracy, the replacement C {as) 1 + as g +as (djv ~ Pn) is allowed, including also 
the remainder terms. The QCD distribution and the shape function are related by a convolution in momentum 
space: 

f (z) = / f dz'dz"S {z - z'z") C {z') f {z") . (73) 
Jo Jo 

The function S]„ reads: 

log E„ = (r - n) log f„ + 9 {n ~ r) [log f „ - log /„] . (74) 
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The explicit expression for n > r is: 



where 



5^ -) 

92^ (w; t) 



(75) 



A2P0 - AiP 



[2 (1 - w - t) log (1 - w - r) - 2 (1 - w) log (1 - w) - (1 - 2t) log (1 - 2t) 



i [2 log (1 - w - r) - 2 log(l - w) - log (1 - 2t)] 



[21og2 (1 - w - r) - 21og2(l - w) - log2(l - 2t)] + 



+^log(l-2r) + §llog(l-w). 



(76) 



\ 

0.8- \ 



Figure 4: Plot of the coefficient function I]„ for the choice of the parameters discussed previously (first 140 
moments). Solid line: NLO distribution; dashed line: LO distribution; dotted line: expansion to order Og of 
the exponent in the NLO distribution. 



Equation (|76D is our main result. Let us comment on it. The terms most affected by long-distance effects are 
those containing the combination 1 — 2w and they are the same in QCD and in the ET. As a consequence, they 
cancel in the difference, proving that the coefficient function is short-distance-dominated to NLO accuracy; we 
believe that this cancellation occurs to any order in PT. The coefficient function is meaningful up to 

1 



n < "crit = ^ exp 



/3o«s (Q2) 



A2 



Note that n'^rit ^i^^ between the critical values for the soft terms and the hard collinear one: 

/ 2 

The coefficient function contains also the jet factor, entirely coming from the QCD form factor: 



Jri 



exp 



/3o 



log(l - w) 



(77) 



(78) 



(79) 
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\ 



Figure 5: Comparison of the coefRcient function with the QCD form factor. Sohd hne: NLO coefficient function; 
big dashed hne: NLO QCD distribution; smah dashed hne: LO coefficient function; dotted hne: LO QCD 
distribution. 

The coefficient function is plotted in fig. 4 for the values of the parameters specified previously. It becomes 
singular at n^rit — 297, i.e. an order of magnitude above the QCD form factor or the shape function. For 
n < 50, the NLO curve has a shape similar to the LO one and it is shifted up with respect to the latter 
by ~ 0.2. In fig. 5 we compare the coefficient function with the QCD form factor. Note that the large- 
approximation produces a cusp singularity in the NLO coefficient function at n = r. 

The computation of the shape function from first principles — namely lattice QCD — requires a large 
amount of theoretical work, consisting of the following steps: 

1. analytic continuation from Euclidean to Minkowski space of the relevant correlation function, a 4-point 
function; 

2. lattice regularization of the effective theory in which the shape function is defined, which involves the 
HQET (Wilson lines off the light cone) and the LEFT (Wilson lines on the light-cone); 

3. computation of the shape function in lattice perturbation theory to full order as; 

4. matching of the lattice shape function with the QCD distribution to full order as; 

5. Numerical computation with a Monte Carlo program of the correlation function on the lattice. 



The first step was taken in ref. |19(| , showing the possibility of a lattice computation of the shape function. Step 
2 has only partially been completed: the lattice HQET has been formulated in ref. |2^ while the lattice LEFT 
— as far as we know — has not been constructed yet. The main difference is that, in the latter case, not 
only soft singularities but also light-cone singularities appear in the euclidean correlation functions; we do not 
attempt at the construction of the lattice LEFT here. No theoretical work exists — as far as we know — on 
points 3, 4 and 5. In general, a lot of work is yet to be done: the purpose of this note was simply showing how 
to resum to all orders the threshold logarithms appearing in the coefficient function. 

Our result (|7^) applies to the coefficient function of the shape function defined in lattice regularization , 
after the identification is done: 

- = c/i, (80) 
a 
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where c is a constant or order 1 and a is the lattice spacing, c is determined, for example, imposing the equality 
of the coefficients of the single logarithms in the lattice shape function and in eq. (steps 3 and 4). In our 
work, we have assumed the consistency of the lattice LEET and the knowledge of the one- loop lattice amplitude. 

In the absence of a (non-perturbative) computation of the shape function, we can compare our coefficient 
function evaluated with a cut-off of a few times the hadronic scale (in practice 1-^2 GeV) directly with the 
experimental data. The mismatch is the non-perturbative component up to a scale ^, namely the shape function. 
The factorization scale fi acts effectively in the coefficient function as a prescription for the Landau pole. In this 
respect, our coefficient function can be considered as a complete QCD form factor, with a prescription to remove 
the Landau pole outside region (Pt|). As we have seen, the singularity is not completely removed, but it is shifted 
to much larger values of N. We stress that our splitting of the QCD form factor into a coefficient function and 
a shape function represents a consistent separation of perturbative from non-perturbative effects; it somehow 
differs from the usual strategy of finding a good prescription for the running coupling in the non-perturbative 
region. 



5 Conclusions 



We have presented a NLO evaluation of the resummed coefficient function for the shape function. Our result, 
together with a one-loop lattice computation that is still missing, allows relating the shape function computed 



with lattice QCD to distributions in heavy- flavour decays in the semi- inclusive region (11). The NLO analysis 



corroborates the results of the leading order one, but it does not reveal any new qualitative feature. The 



coefficient function is short-distance-dominated in the relevant region (11), and corrections to factorization are 
expected to be of the order A//i, i.e. to involve one inverse power of the factorization scale. The NLO coefficient 
function contains also a jet factor that takes into account the emission of hard coUinear gluons. The latter 
process is indeed not described by the effective theory, which takes into account only soft emission up to the 
scale fi. On the quantitative side, NLO effects increase the coefficient function by ~ 20% for /z = 2 GeV and 
n < 50 in the case of B decays. 

In general, the process (0) contains three different scales:^ 



i) » A^ 

a) Qk+ - QA > A^, 
in) ~ A^. (81) 



The emission of gluons with transverse momenta ranging between i) and ii) is reliably computed in perturbation 
theory, while fluctuations in region Hi) ask for a non-perturbative treatment. With factorization by means of 
the shape function, another scale is introduced in the problem, intermediate between ii) and Hi) : 

iv) ^k^ ^ fi A > A^. (82) 

The fluctuations between scales Hi) and iv) are factorized in the shape function. The general idea behind 
the shape-function approach is that soft contributions are more "non-perturbative" than coUinear ones in 
double- logarithmic problems in region (0), because they involve smaller transverse momenta. These soft non- 
perturbative effects can be represented as hadronic matrix elements of non-local operators involving Wilson 
lines. For a given process, one identifies a dominant kinematical configuration and replaces the hard partons 
with Wilson lines with the same momenta. 



To summarize, we have presented a NLO evaluation of the resummed coefficient function of the shape 
function, which confirms the consistency of the effective theory approach and allows a more accurate theoretical 
analysis. 

wish to thank G. Korchemsky for a discussion on this point. 
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